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We investigate the behaviour of tensor fluctuations in Loop Quantum Cosmology, focusing on a 
class of scaling solutions which admit a near scale-invariant scalar field power spectrum. We obtain 
the spectral index of the gravitational field perturbations, and find a strong blue tilt in the power 
spectrum with nt ~ 2. The amplitude of tensor modes are, therefore, suppressed by many orders of 
magnitude on large scales compared to those predicted by the standard inflationary scenario where 
nt ~ 0. 

PACS numbers: 



I. INTRODUCTION 

A key result of the inflationary scenario of the very 
early universe [l| is that it gives rise to a stochastic back- 
ground of gravitational wave radiation (tensor perturba- 
tions) 0. Such a background is in principle observable 
and could be used to distinguish between different models 
of inflation [3|. Alternative proposes such as the ekpy- 
rotic/ cyclic scenario 0, Hj or phantom super- inflation 
[1, Q lead to very different predictions for the spectral 
tilt of tensor perturbations when compared with the sim- 
plest inflationary models 0, implying the gravitational 
wave background could also be a powerful discrimina- 
tor between competing theories. In this work we will be 
concerned with the gravitational wave background pro- 
duced by super-inflationary scenarios within Loop Quan- 
tum Cosmology (LQC). 

Loop Quantum Gravity (LQG) is a background in- 
dependent and non-perturbative canonical quantisation 
of general relativity based on Ashtekar variables: su(2) 
valued connections and conjugate triads. The variables 
used in the quantisation scheme are holonomies of the 
connection and fluxes of the triad. Restriction to sym- 
metric states followed by quantisation using the tech- 
niques of LQG gives rise to LQC Although much 
regarding LQC is not fully understood, in particular its 
relation to full LQG, it has produced a number of intrigu- 
ing results and resolved the problem of singular evolution 
present in the earlier Wheeler de Witt approach to quan- 
tum cosmology [l0| . 

One particularly useful approach to studying the early 
universe in the context of LQC is that of deriving and 
studying effective equations of motion [Til [T^ [isl. [isj . 
In the context of inflation in LQC many in triguing ef- 
fects have been investigated in this setting [l3l flGlflTl. 
[isj . While such equations cannot probe fully quantum 
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regimes, they nevertheless incorporate quantum modifi- 
cations into classical evolution equations whilst avoiding 
the interpretational difficulties inherent in fully quantum 
equations. In isotropic settings two sets of modifications 
have been predominately considered (see Ref. [l^ for 
a discussion of other corrections which may also need 
to be included in certain regimes). The first originates 
from the spectra of quantum operators related to the 
inverse triad [ll|, while the second arises from the 
use of holonomies as a basic variable in the quantisa- 
tion scheme [l^. One point to note here is that in the 
regime of very small scale factor, the discrete structure 
of space is so strong that any inhomogeneous conflgura- 
tion would be far from being isotropic. Therefore, despite 
being able to write the effective perturbation equations 
for a purely isotropic scenario in this regime, one should 
keep in mind that in reality care needs to be taken re- 
garding the cosmological interpretation of perturbations 
and their evolution [2l|, [g^l • 

The two modifications have rather different origins, yet 
both can give rise to a period of super-inffation during 
which the Hubble rate rapidly increases, rather than re- 
maining nearly constant as is the case during standard 
slow-roll inffation. Since the relative status of the two 
modifications discussed is at present unclear, most works 
have taken the pragmatic approach of studying the dy- 
namics when each of the modifications is considered sep- 
arately, but not including both sets of modifications si- 
multaneously. 

In earlier work [TtI . [TsI , we explored super- inflationary 
scalar field models in the presence of each set of modifi- 
cations in turn. We worked in an approximation which 
neglected back-reaction, and by studying scaling, power- 
law solutions we derived the form of the potential re- 
quired in each case for perturbations to the scalar field, 
produced during the super-inflationary phase, to have 
a scale invariant spectrum. Moreover, while the period 
of super-inflation in each case must be short lived, we 
showed that this is not a barrier for these scenarios since 
only a small number of e-folds of super-inflationary evo- 
lution is required to solve the horizon problem. In that 
work we did not, however, consider the spectrum of ten- 
sor perturbations produced by these scenarios, and this 
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is the question to which we now turn. This question is 
particularly timely since the effect of the two kinds of 
modifications on the evolution of tensorperturbations in 
LQC has only recently been derived [23]. The modifi- 
cations have already been utilised to calculate both cor- 
rections to the gravitational wave spectrum if standard 
inflation occurred in the presence of LQC corrections [23| 
and the gravitational wave spectra from the 'big bounce' 
[25| which can occur in LQC. 

Similar calculations have already been carried out and 
reported in the literature for different types of cosmolo- 
gies. In particular, gravitational wave perturbations 
for ekpyrotic models of a collapsing universe 0, and 
for phantom superinflation scenarios have been com- 
puted. Despite being based on different physical concepts 
and behaviours, both cosmologies predict a strongly blue 
tilted spectrum of tensor perturbations, and since these 
have not yet been observed, it suggests that in both sce- 
narios, the amplitude of these fluctuations are suppressed 
on large scales by many orders of magnitude compared 
to those predicted by standard inflation. This is not un- 
expected, as it has been pointed out in Ref. [2^ that a 
duality exists between the ekpyrotic collapse and the dy- 
namics of a universe sourced by a phantom field. More- 
over, a scale factor duality maps the ekpyrotic collapse 
onto the superinflationary scaling solutions in LQC [27|. 

Given the connections found between LQC and the 
ekpyrotic and phantom scenarios, one might expect that 
LQC also predicts a large blue tilted spectrum for the 
tensor perturbations. However, as the connections are 
at the level of the background equations and since LQC 
corrections also arise in the evolution equations of tensor 
perturbations themselves, this expectation may not be 
realised and a careful analysis is necessary to confirm it. 
We will see below that these corrections do not spoil this 
conjecture after all. 

The structure of the paper is as follows. In Section 
im we discuss the inverse triad modifications. First we 
review the background dynamics which give rise to a scale 
invariant spectrum of scalar field perturbations, and then 
we consider the evolution of tensor perturbations in this 
setting calculating their spectrum. We then repeat the 
exercise for holonomy corrections in section IIIII Finally 
we conclude in section HVl 
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state their basic properties. We note that we are implic- 
itly considering either positively curved or topologically 
compact models. This ensures that the size of the fidu- 
cial cell does not enter in the equations of motion. D and 
S are both functions of the scale factor, and their form 
changes depending on the values of two ambiguity param- 
eters: I which takes values in the range < I < 1, and 
j which takes half integer values. When the scale factor 
approaches zero, D and S also approach zero, whereas as 
a increases above the critical value a*, which depends on 
J, they both tend to unity. As noted in the introduction, 
for small values of the scale factor a <C a*, the assump- 
tion of isotropy may be violated. Whilst recognising the 
careful attention this regime deserves, and acknowledg- 
ing more work is required to clarify the validity of the dy- 
namical equations we use together with the assumption 
of isotropy, we hope to pave the way by demonstrating a 
method of calculation which can easily be employed once 
new light is shed on currently uncertain sections of the 
theory. 

The unperturbed isotropic equations of motion take 
the form of a modified Friedmann equation 
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where a dash represents differentiation with respect to 
conformal time and we have assumed that any curvature 
contribution has become subdominant, and the scalar 
field equation is 
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TENSOR DYNAMICS WITH INVERSE 
TRIAD CORRECTIONS 



For convenience, from here on we will work in units in 
which = 8nG = 1. 

In this study our interest is in the evolution of ten- 
sor perturbations about this isotropic background. Ten- 
sor perturbations are defined as the transverse and trace 
free part of the perturbed spatial metric, and represent 
gravitational wave perturbations propagating on the un- 
perturbed background spacetime. They can be further 
decomposed into two polarisation modes represented by 
X and -|-, and in LQC, with inverse triad modifications 
included, the equation of motion for these modes has re- 
cently been derived to be [1^ 



We first consider the cosmological equations which fol- 
low from including modifications associated with the in- 
verse triad in LQC. This introduces two functions into 
the dynamics, Dji{a) and Sj{a). These functions arise 
because of the presence of powers of the inverse scale 
factor in the Hamiltonian constraint for an isotropic and 
homogeneous universe. A full discussion of the origin of 
these terms can be found in Ref. [Tp . [T^ (a summary can 
be found in appendix B of Ref. [28[), but here we simply 
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h'^^- S^W^hy,^+=Q . (3) 



hx.+ is a tensorial quantity, but from here on to avoid 
clutter we will drop the x and -I- subscripts, and take 
h to represent the magnitude of one of the polarisation 
modes, but always keeping in mind that both modes are 
present. 
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The background power law solution and scale 
invariant scalar field dynamics 



B. The primordial spectrum of tensor 
perturbations 



In Refs. [I^ [27[, it was shown that in the regime 
a <C a* there exists a power law solution to the equa- 
tions of motion (IT])-([2]) which is a stable attractor to 
the dynamics. In this regime D{a) ~ D^a^, with 
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takes values in the range 9 < rt < oo. The function S{a) 
may be similarly approximated by S{a) ~ S'*a'", where 
Si, = (3/2)a^'" and r = 3, though we keep r arbitrary in 
our calculations for generality. For a ^ a* , S** « « 1 
and r = n = 0. The power law solution for a ^ a* exists 
for negative power law potentials of the form V = Vo(j)^ 
and gives rise to the dynamics 



a(r) = Ai-Tf, 
cp'ir) = 
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where for an expanding universe t is negative and in- 
creasing towards zero. A is an arbitrary normalisation 
constant, a = l — n/6 with /3 andp being related through 
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One can expand about this solution in terms of fast roll 
parameters, in order to generalise the potentials which 
can be considered [l8l |. but here for simplicity we will 
consider only this exact solution. 

For — 1 < p < the solution given above represents a 
universe undergoing super-inflationary evolution during 
which H = a/a {a, dot denotes differentiation with re- 
spect to cosmic time), the Hubble rate, increases. A par- 
ticularly interesting case then occurs as p tends to zero 
from below. This represents a universe in which the scale 
factor is almost constant, but H increases rapidly. Con- 
sidering scalar field perturbations about the background 
field, it was shown in Refs. [13) EB], that the spectrum of 
scalar field perturbations attains scale invariance in this 
limit. Moreover, because H increases so rapidly, the hori- 
zon problem is solved during this phase with only a small 
number of e-folds required. This raises the intriguing pos- 
sibility that these perturbations could be responsible for 
the observed CMB anisotropics and hence for structure in 
the universe. If this were the case, no period of standard 
inflation in which H remains nearly constant for roughly 
60 e-folds of expansion would be required (where num- 
ber of e-folds is defined as TV = ln(a/ai)). A natural, and 
indeed important question is: 'what is the spectrum of 
primordial tensor perturbations which accompanies this 
scale invariant spectrum of scalar field perturbations?'. 
This is the question to which we now turn. 



To calculate the spectrum of gravitational waves pro- 
duced during super-inflation we follow the standard pro- 
cedure. Noting that 2a'^h is canonically conjugate to h' / S 
such that m 
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the system is quantised by promoting h and h' to op- 
erators and the Poisson brackets to commutators. We 
have 
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h is then decomposed into Fourier modes 



(27r)3/2 



hkiT)ake''^-'' + hl{T)ale 



k 



(9) 



(10) 



where, considering Eq. ([3]), we see that each mode hk 
obeys the evolution equation 
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The power spectrum for one polarisation state of ten- 
sorial fluctuations is given by the standard expression 
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In order to evaluate Eq. (fT^ we must solve Eq. pip . 
Considering the power law solution for the regime a <C a^, 
(Eq. (j4|)), and using the form of S in this regime, we find 
that Eq. (fTT|) becomes 
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and we have normalised the solution such that at and aj^. 
satisfy the usual raising and lowering operator algebra 
while h and h' satisfy the commutation algebra such 
that only the forward moving solution is selected in the 
asymptotic past (the adiabatic vacuum). The solution 
p4)) has the expected behaviour that each k mode be- 
gins in an oscillatory state where normalisation occurs, 
and evolves into a non-oscillatory state once each given k 
mode crosses a suitably defined horizon. From the form 
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of ((T4)) . it is clear that for tensor modes horizon cross- 
ing occurs when k « —aH/pS. A given mode can only 
be considered to become a classical perturbation once it 
crosses this horizon. 

Taking the solution ([14]) and employing (fT2|) we find 
that 



22(.-l) p(^)2 

^2 r(3/2)2 



(16) 



We have evaluated ([H]) in the limit where the modes 
are outside the defined horizon (i.e. x — > 0), and have 
accounted for both polarisation states. The mode ko that 
corresponds to the largest scales on the CMB and fee is 
defined as the last mode to cross the horizon at the end 
of super-infiation (fco = Ho/Sc ~ He)- 

A couple of observations are in order. The first is that 
in the limit of interest, p — > 0, the spectrum is blue tilted 
with a tensor spectral index given hy nt —2, where nt 
is defined, as usual, by Ph oc fc"' . This implies that on 
large scales the spectrum is hugely suppressed. The sec- 
ond important point is that the magnitude of the spec- 
trum is fixed by the Hubble rate at the end of inflation. 
This also fixes the magnitude of scalar perturbations and 
ultimately has to be normalised such that the scalar per- 
turbations have the correct magnitude to account for the 
CMB anisotropics. In the scenario at hand such a nor- 
malisation is difficult to determine since the scalar field 
perturbations (which are close to scale invariant) must 
be related to curvature perturbations (as discussed at 
length in Refs. [3), and this step cannot be performed 
at present. Nevertheless in the following section we will 
make the reasonable assumption that He « 10~^ cor- 
responding to the GUT scale in order to calculate the 
present-day spectrum of gravitational waves produced by 
this super-inflationary scenario. Furthermore, we will see 
that the conclusion - that the spectrum is unobservably 
small - is insensitive to the choice of normalisation within 
reasonable bounds. 



C. The present-day spectrum 

At the end of super-inflation all classical perturbations 
modes are outside the horizon. We will assume that re- 
heating occurs instantaneously at the energy scale He, 
and hence that the universe becomes radiation domi- 
nated at this point. Moreover we will assume the universe 
is classical after reheating, and any quantum corrections 
(similar to D or S) are absent from the dynamics. From 
this point onwards modes will begin to re-enter the cos- 
mological horizon. We schematically illustrate this dy- 
namics in Fig. [T] 

To convert the primordial spectrum, Eq. to the 

spectrum which would be observed today we employ the 
numerically obtained transfer function [2^ [30| 
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FIG. 1: Schematic illustration of the evolution of modes k 
which exit the horizon (thick solid line) during super-inflation 
and re-enter during the standard radiation or matter era. k^^ 
is the scale corresponding to the size of the observable uni- 
verse. 



where fcg = HqQo, is again the k mode corresponding 
to the largest scales today, and k^q — a^qHoq, where 
eq stands for quantities at radiation matter equality. 
This form arises from the observation that tensor modes 
outside the horizon are roughly time independent while 
modes inside decay as a~^, together with the evolution 
of H (X a~2 during radiation domination and H oc a^'^/^ 
during matter domination. 

Using this transfer function we can calculate the 
present day spectrum of tensor perturbations once we 
flx He and kc/k^. A sensible estimate for He is the 
GUT scale, while an absolute upper limit is given by 
the Planck scale. If we also make the assumption that 
as many modes exit the tensor horizon as scalar modes 
exit the scalar horizon, k^/k^ can be fixed by the require- 
ment that the horizon problem be solved. In Ref. [l^ 
we found that this required ko/k^ w e~®°. 

A useful physical quantity that can be used to express 
the present day spectrum of gravitational waves is ^gw, 
the gravitational wave energy per unit logarithmic wave 
number in units of the critical density, pcrit 
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Figure [2] shows a plot of the present tensor abundance 
for a number of choices of the parameter p, together with 
the strongest observational constraints from current and 
future experiments searching for gravitational waves [sH . 

It can be verified that even the limiting value of H^ ~ 1 
(giving the largest spectrum possible) , leads to an unob- 
servably small spectrum of tensor perturbations. 
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FIG. 2: Predicted present abundance of gravitational waves 
produced during the super-infiationary phase of LQC with in- 
verse volume corrections assuming that = 10~^ in Planck 
units. The solid line corresponds to the standard inflation 
abundance (n = r = 0). Also indicated are the present 
bounds and future sensitivities. 



III. TENSOR DYNAMICS WITH HOLONOMY 
CORRECTIONS 

We now turn our attention to the second set of effec- 
tive equations, those which arise from considering that 
holonomies are the basic operators to be quantised in 
LQC. 

The isotropic unperturbed dynamics is described by 
the Friedmann equation p^ : 



1 2 
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which is modified from the classical equation by the in- 
clusion of a —p^ term, where cr = 3/ (27^) is a constant in 
an exactly isotropic model, and where 7 is the Barbero- 
Immirzi parameter [s^ [3^ . Once again we are assuming 
either a flat universe or that the curvature contribution 
can be safely neglected. Our interest here is in a scalar 
field dominated universe, and hence p = (j)'"^ /2a? + V{(t)). 
The scalar field equation of motion remains unchanged 
from its classical form 



+ a^V.^ = . 
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We stress that as we are studying inverse volume and 
quadratic corrections separately, we do not include the 
D and S functions in the equations of motion. 

The form of the evolution equation for tensor pertur- 
bations when holonomy corrections are included has re- 
cently been derived to be [1^ 
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are quantum corrections to the classical dynamics that 
become unimportant when p <C 2a. From here on we 
again drop the x and -I- subscripts as we did in the inverse 
triad case. 



A. Power- law solution and scale invariant scalar 
field perturbations 

We are interested in high density regimes where p ap- 
proaches the bounding value of 2a. In this case, the 
term within brackets of Eq. tends to zero and the 
behaviour of the equations alters significantly compared 
with the classical behaviour. Indeed in this regime we 
have H > and for an expanding universe super-inflation 
takes place. In our previous work p^ . we showed that in 
this regime there exists an approximate power-law solu- 
tion for a potential of the form V = 2a — UQe"^"^ . More- 
over, this solution is an attractor as we demonstrated 
both analytically and numerically [11] . The full solution 
is given by 



a(T) = A(-t)p, 



h'{T) = V-2p(P+l) 



VIt) = 2a- 



P(2p-1) 
{arf 



(24) 
(25) 

(26) 



where A is an arbitrary normalisation constant. A and p 
are related through 



X' = ~2 
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Once again — 1 < p < corresponds to a universe un- 
dergoing super-inflation, and moreover the limit p ^ 
from below leads to a scale invariant spectrum of scalar 
field perturbations. Furthermore, as was the case for the 
super-inflationary solution we studied in the presence of 
inverse volume corrections, only a small number of e-folds 
are required to solve the horizon problem. In our previ- 
ous work we generalised the form of the potential that 
could be considered by expanding about this solution, 
but for simplicity in this work we will consider only this 
exact form. 

We now turn our attention to the question of what 
spectrum of tensor perturbations accompanies the scale 
invariant spectrum of scalar field perturbation in this ver- 
sion of super-inflation. 
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B. The primordial spectrum of tensor 
perturbations 

We again follow the standard procedure for calculating 
the spectrum of tensor perturbations. In this case a^h is 
canonically conjugate to h'. To quantise the system h 
and h' are promoted to operators and Fourier decom- 
posed according to Eq. PH)) . Each /i^. mode satisfies the 
evolution equation 



h'l + 2-h' 



,,2 Mp + 1) 



during the scaling solution, where we have used p 
and Eq. 

The solution to Eq. is given by 



-^/^ff(l)(-fcT), 



where 
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where we have normalised the solution by the require- 
ment that Ok and aj^ satisfy the usual raising and lowering 
operator algebra while a^h and h' satisfy their commuta- 
tion algebra, with only the forward moving solution being 
selected in the asymptotic past (selecting the adiabatic 
vacuum). Finally, utilising Eq. and evaluating this 
expression in the limit that the modes are outside the 
horizon (fc < —Ti./p in this case), leads us to the same 
expression for the primordial power spectrum (|16p . The 
spectral index can clearly be seen to be given by nt = 2 in 
the limit p ^ 0, and the amplitude fixed by He. We note 
here that in cases where the dynamics of tensor pertur- 
bations could be represented by the standard equation of 
motion and the evolution of the scale factor is given by 
Eq. l|24|) . having p ^ will inevitably result in rit = 2. 
This has indeed been shown to be true explicitly for the 
collapsing ekpyrotic scenario Q and for the evolution of 
a universe sourced by a phantom field Q- At this point, 
making some reasonable assumptions, we can proceed, 
as we did for the solution in the inverse triad case, to 
calculate the present day spectrum of gravitational wave 
perturbations using the transfer function Eq. (|17p . We 
show the abundance of tensors in Fig. [3] and note that 
the result is very similar to the inverse volume case. 



IV. DISCUSSION 

We conclude that super-inflation in both versions of 
the corrections in LQC predicts a strong blue spectrum 
of gravitational waves, hence, their abundance is strongly 
suppressed on the large scales and it is many orders of 
magnitude smaller than the value predicted in the stan- 
dard inflationary scenario. 
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FIG. 3: Predicted present abundance of gravitational waves 
produced during the super-inflationary phase of LQC with 
holonomy corrections. The solid line corresponds to the stan- 
dard inflation abundance. Also indicated are the present 
bounds and future sensitivities. 



It is now important to discuss our results in the light of 
other investigations in the literature. In particular, it is 
claimed in Ref. [s^l that the abundance of gravitational 
waves generated during super-inflation under inverse vol- 
ume corrections is above the current bounds. We note 
however that the authors did not consider the evolution 
of the background in a scaling solution and the S(a) cor- 
rection was not used. Moreover, the full expression for 
the tensor perturbations was obtained more recently [1^ 
and therefore, it was not used in that work. In a more 
recent analysis [IJl focusing on the holonomy corrections, 
it is also found, like in our work, that the spectrum of 
gravitational waves must be blue. However, a scaling so- 
lution was not used and the expansion of the universe 
was assumed to be close to de Sitter, hence, a direct 
comparison with our work is in fact not possible. 

One needs to be cautious about making general con- 
clusions in the context of LQC, as the theory is currently 
far from complete. As mentioned previously, the intro- 
duction of inhomogeneities in the small a regime is likely 
to break the assumption of an isotropic universe. At- 
tempts are being made to gain better understanding of 
this region [2l[. Moreover, there is the possibility that 
higher order perturbative correction terms could play a 
role, or that quantum backreaction might significantly 
modify the background dynamics [l^ . 

Our work on super-inflationary scenarios in LQC also 
has a number of other possible drawbacks. We have 
treated the inverse triad and holonomy corrections sepa- 
rately while they should be dealt with together in a real- 
istic set up. Though the fact that both sets of corrections 
lead to such similar phenomenology gives us some reas- 
surance that combining them could lead to qualitatively 
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similar results. Furthermore, so far we have not investi- 
gated the evolution of scalar metric perturbations. The 
derivation of the full equations for these is still in progress 
[35| and although these equations are not required for the 
calculation of tensor perturbations presented here, they 
are required to understand how the scale-invariant scalar 
field perturbation is related to the observed curvature 
perturbation. Finally, we should mention recent work 
where it has been shown that the behaviour of the LQC 
equations with holonomy corrections in the presence of 
negative exponential scalar field potentials leads to sud- 
den singularities where the Hubble rate is bounded, but 
the Ricci curvature scalar diverges [sl]. Given that for 
the case of holonomy corrections our superinflation sce- 
nario requires a scalar field potential with a negative ex- 
ponential part, this serious problem needs to be avoided. 
In the scenarios we consider, however, the potentials only 
need to be of the form which gives rise to the power-law 
behaviour while superinflation is taking place. After this 
phase of evolution the potential can change in form. For 



example, any potential which tended to zero after the 
field evolved through the region which gives rise to the 
super-inflation phase would avoid the sudden singularity 
problem. 

Despite the ambiguities and uncertainties in the the- 
ory, if we were to make observational claims for the scal- 
ing solutions we have considered in the current state of 
LQC, it would be that if gravitational waves are observed, 
they would rule out this scenario of superinflation in LQC 
sas it stands. 
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